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We study the structure characteristics of complex networks from the representative eigenvectors of 
the adjacent matrix. The probability distribution function of the components of the representative 
eigenvectors are proposed to describe the localizations on networks where the Euclidean distance 
is invalid. Several quantities are used to describe the localization properties of the representative 
states, as the participation ratio, the structural entropy, the probability distribution function of the 
nearest neighbor level spacings for spectra of complex networks. The whole cellular networks in 
real world, the Watts-Strogatz small world and Barabasi- Albert scale-free networks are considered. 
The networks have nontrivial localization properties due to the nontrivial topological structures. It 
is found that the ascend-ranked series of the occurring probabilities at the nodes behave generally 
multi-fractal. This characteristic can be used as a structure measure of complex networks. 
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I. INTRODUCTION 



Recent years have witness an avalanche of research on 
complex network and its applications in diverse fields 
[l[ . Structure measures of complex networks are the cor- 
nerstone to understand relations between the structures, 
dynamics and functions. Real world networks generally 
have nontrivial properties as the small- world |2jj, scale- 
free 0, motif jj], modularity, hierarchy Q, fractal @ 
and so on. The small-world effect is that in average the 
nodes can reach each other with only a small number of 
hops. The scale-free refers to the number of edges per 
node obeys a right-skewed distribution. It is also found 
that some special subgraphs containing several connected 
nodes, called motifs, occur with significant probabilities 
compared with that in the corresponding randomized 
networks. These three individual, pair or local pattern- 
based properties are called micro-properties. On the 
other hand, the modularity is a kind of macro-property 
represents that a network can be separated into loosely 
connected groups within which the nodes are densely con- 
nected, respectively. 

To a certain degree, dynamics on networks can be re- 
garded as the transport processes of mass, energy, signal 
and/or information at different structure scales pi [8|. 
Sometimes we have to deal with networks with unrea- 
sonable large number of nodes and edges, e.g., the neu- 
ron networks and the World- Wide- Web networks, when 
designing a coarse-grain procedure is essential [9|. The 
patterns at different scales may provide a reasonable so- 
lution to these problems. It is found that some real world 
networks have hierarchical structures, in which the small- 



world and the scale- free properties can coexist @ . More- 
over, many real world networks behave self-similar at dif- 
ferent structure levels (fractal) @ . 

Though great progresses have been archived, the mea- 
sures of complex networks are not yet fully understood. 
Just as pointed out by Newman [Toj , that our techniques 
for analyzing networks are at present time no more than 
a grab-bag of miscellaneous and largely unrelated tools, 
and we still do not have a systematic program for charac- 
terizing network structures. Furthermore, the measures 
of network structures, such as the micro-properties, the 
patterns at different scales and the macro-properties, are 
generally a simple application of the concepts in graph 
theory, bioinformatics, social science and fractal theory, 
namely, they are not dynamics-based. We can not ex- 
pect simple and reasonable relations between the mea- 
sures and the dynamical processes on networks. 

The lack of powerful tools to characterize network 
structures is an essential bottleneck to understand dy- 
namical processes on networks. One typical example is 
the synchronizabilities of complex networks. Detailed 
works show that almost all the structure measures affect 
the synchronizabilities in complicated ways [ll|, based 
upon which we can not reach a clear picture of the mech- 
anisms for synchronization processes on networks. 

Dynamics-based measures of complex networks may 
be the key to the problems. The structures of complex 
networks can induce nontrivial properties to the physical 
processes occurring on them. The physical processes in 
turn can be used as the probe to capture the structure 
properties. Well studied dynamical processes, such as the 
random walks [H|,[l3| and the Boolean dynamics [l4j], can 
be good candidates as probes. To cite an example, the 
random walks on complex networks that biased towards a 
target node show a localization-delocalization transition 
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In the present paper, we map networks to large clus- 
ters, namely, the nodes and edges to atoms and bonds 
between them, respectively. The localization properties 
of electrons in the clusters can be used as measures of the 
structure properties of the networks. We try to detect 
the global symmetries from the spectra and the eigen- 
vectors of complex networks. Very recently, much atten- 
tions have been focused on detecting global character- 
istics embedded in spectra of complex networks due to 
their potential application in understanding the organi- 
zation mechanisms and the synchronization dynamics of 
complex networks @, [l5|, [llj] . To our best knowledge, it is 
the first time to detect the global characteristics of com- 
plex networks from the eigenvectors which contain more 
information about the system than eigenvalues. 

Besides as a measure of network structures, the 
structure-induced localization may have potential appli- 
cation in understanding the electronic properties of ma- 
terials such as conductive polymers and carbon nanonets. 
The intrachain windings in conductive polymers can in- 
troduce long-range edges into the original one dimen- 
sional (1-D) systems, resulting in nontrivial network 
structures [17|, LL8J- ^ is also found that random net- 
works of carbon nanotubes, called nanonets, can mimic 
a variety of basic electronic functions from the conductive 
properties of metals to the less conductive characteristics 
of semiconductors (19J. Indeed, nanonets have paved the 
way for the carbon to serve as the foundation for future 
electronic devices. The effect of network structures on 
electronic properties is one of the most active topics in 
recent years 0, [2(| . 

The paper is organized as follows: In Sec. II the concept 
of localization on complex networks is presented. The 
occurring probabilities on the nodes are proposed to de- 
scribe quantitatively the localization effects. In Sec. Ill 
the methods to measure the localization properties are 
described in detail. The participation ratio, the struc- 
tural entropy and the probability distribution function 
of the nearest neighbor level spacings of spectra are used 
to illustrate the localization in a global way. The wavelet 
transform is then used to find the detailed structure prop- 
erties of the probability distribution function of the oc- 
curring probabilities on the nodes. As examples, we con- 
sider the Watts Strogatz small-world and the Barabasi- 
Albcrt scale-free modeling networks and the whole cellu- 
lar networks in real world. The results are shown in Sec. 
IV. We will be shown that the global symmetries in net- 
works can induce multi-fractal structures in the eigenvec- 
tors. As a conclusion, the nontrivial structures of com- 
plex networks can induce significant localizations, which 
in turn can be used as a global measure of the structure 
symmetries. 



II. THE LOCALIZATIONS ON NETWORKS 

We consider a undirected complex network with N 
identical nodes, whose topological structure can be de- 



scribed by an adjacent matrix A. The elements Aij are 
1/0 if the nodes i and j are connected/disconnected, re- 
spectively. If we consider the nodes as atoms and the 
edges as bonds, the network can be mapped to a large 
molecule [lq j. For an electron moving in such a molecule, 
the tight-binding Hamiltonian is, 

JV N 

H = ^e n - \n) (n\ + ^ A mn • t mn ■ \m) (n\, (1) 

n— 1 m^n 

where e n is the site energy and t mn ■ A mn the hoping 
integral for the bond between sites m and n. 

A tight-binding Hamiltonian of Eq.(l) is usually used 
to study the disorder-induced localizations. In the 
present form, the matrix A is explicitly introduced to de- 
scribe the structure of the system. For a one dimensional 
(1-D) perfect regular lattice, we have s n — s,t rnn = t 
and A mn = 5(m — nil). The Bloch wave function of an 
electron extends all over this perfect regular lattice. Dis- 
order structures can induce a transition from extended to 
localized states. The wave function for a localized state 
decreases exponentially with the distance from its center. 
The disorder effects include the random distributions of 
the site energies (e„), the hoping integrals (t mn ) and the 
edges in structures (A mn ). The disorders come from the 
different kinds of atoms on the lattice points, the differ- 
ences of the separations of successive lattice points and 
the randomness in structures. At the same time, there 
may be some symmetries in the distributions of the site 
energies, the hoping integrals and the edges, which may 
lead to derealization of the wave functions. 

In the usual Anderson model [2l|, the disorder ef- 
fect due to the random distribution of the site ener- 
gies is considered, i.e., e n is a random variable satis- 
fying a certain distribution probability function while 
tmn = t, A rnn = 8(m — n ± 1). The site energies may obey 
a special distribution rather than that in the Anderson 
model, as a periodic [22j or a power-law [23| function. In 
literatures 0, [H[ , a one-dimensional quasicrystal model 
is introduced that the separation of two successive lat- 
tice points takes one of the two values u and v. This 
model considers the disorder effect of the distribution of 
the hoping integrals. We have e n =const., t mn = t u or 
t v and A mn = 5(m — n ± 1). t u and t v are the hoping 
integrals corresponding to the separations u and v, re- 
spectively. It is found that quantum systems with quasi- 
periodic structures will be in an intermediate state, which 
can be described with critical wave functions. A critical 
wave function obeys a power-law with respect to the dis- 
tance from its center. 

To investigate the problems as vibration spectra of 
glasses , instantaneous normal modes in liquids, electron 
hopping in amorphous semiconductors and combinatorial 
optimization, Euclidean random matrix (ERM) models 
are widely used in literatures (26[, in which the disorder 
is due to the random positions of the sites, and the ma- 
trix elements are given by a deterministic function of the 
distances. 
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The models mentioned above generally focus on the 
disorder effects of the site energies and the hopping inte- 
grals. These models have also been extended to nontriv- 
ial structured systems such as the Cayley tree [27j and 
the small- world networks [17|. Nontrivial effects of the 
structures of the systems are reported, but the interplay 
between the disorders due to the site energies and that 
due to the structures makes it difficult to distinguish the 
structure disorder effect from the site energy disorder ef- 
fect. 

In the networks considered in the current paper, how- 
ever, the nodes are all identical and the disorder effect 
comes from the nontrivial topological structure. We fo- 
cus our attentions on the disorder effect of the network 
structure, that is, we assign e„ = and t mn — 1, which 
leads H = A. The localization on the network refers to 
the network structure- induced characteristics of the wave 
functions for this system. The usual Anderson model [2l[ 
is a special case that in the networks there exist connec- 
tions only between the nearest neighbors in Euclidean 
space. 

Statistically, the structures of networks should display 
certain symmetries due to the general rules obeyed in 
the construction of the networks. Recent works demon- 
strate that many theoretical and real world networks 
have statistically self-similar structures [f|. Therefore, 
there are two competitive mechanisms determining the 
wave function property, the randomness of the bonds in 
the netowkrs tends to cause localization of wavefunction, 
whereas the symmetries of the networks intend to make 
wave function extended. We thus expect rich structures 
embedded in the wave functions. As it is well known, ape- 
riodic crystals lead to the fractal wave functions [25|, [28| ■ 
An interesting question is then, how the global symme- 
tries of networks affect the localization properties. The 
localization can be used as a probe of the characteristics 
of the network structures. 

The states in the center of the energy band have the 
best chance to remain as extended for a moderately dis- 
ordered system. The eigenvector corresponding to the 
special eigenvalue close to the center of the spectrum 
for a network, denoted by E c , is employed as the rep- 
resentative state to illustrate the characteristics of the 
considered system. 

In the traditional study of wave function localization, 
the physical systems have deterministic structures in real 
world Euclidean space, which leads to natural definitions 
of the localized, intermediate and extended states of the 
systems. Obviously, these definitions are invalid for gen- 
eral complex networks without deterministic structures 
in Euclidean space. In this paper, we describe the local- 
ization effects with the probability distribution function 
(PDF) of the occurring probabilities at the nodes, i.e., the 
values of the components for the representative eigenvec- 
tor. Based on the PDF of the occurring probabilities, 
the traditional definitions are extended to a much more 
general version to describe the localization properties on 
complex networks. 



In Euclidean space, for a state \l/ (r), the occurring 
probability is p(r) = \^{r)\ = F(r). Because the value 
of the distance r distributes homogeneously in the con- 
sidered region, we can regard it as a homogenously dis- 
tributed random variable. The direct sampling method 
in Monte Carlo simulations tells us that the probability 



distribution of p should be P{p) oc 



dp 



Hence, it is 



reasonable to define the localized, critical and perfectly 
extended states on complex networks with the PDFs of 
the occurring probabilities, 



S(p-Po), 
P(p) oc <( p-V+n) |ry > , 

~ 77 = 0, 



-,-(1+^) 



extended 
critical 
localized 



(2) 



The PDF of the representative function is a very powerful 
measure to capture the localization properties. It can 
be used to find the localization properties without using 
distance in real world Euclidean space. 

Because no derivative exists for a fractal wave-function 
in Euclidean space, the extension procedure in defining 
critical and localized states on networks can not be sim- 
ply used to define fractal property on networks with the 
PDF of the occurring probabilities. In the present pa- 
per, we detect directly the fractal characteristics in the 
ascend-ranked series of the the occurring probabilities, as 
described in detail in Section 111(C). 



III. METHODS 

A. Structural Entropy 

We denote the representative state with V = 
{Vx, V2, ■ ■ ■ , V/v). The occurring probabilities at the 
nodes are p m = \V m \ ,m = 1,2,-- - ,N. The localiza- 
tion extent of the state can be described quantitatively 
with the participation ratio [2!| [3(| , 



(3) 



For a perfect extended state we have Q = 1, while for 
a state strongly localized on one node it tends to ■h. 
Generally, Q should be in the range of l] . 

However, this participation ratio can capture only the 
primary-level complexity in the localization properties, 
namely, the extension of the representative eigenvector to 
N ■ Q nodes on the network. Many PDFs with different 
localization behaviors may result in the same Q. The 
simplest one is a step-like function that on TV • Q nodes 
the occurring probabilities are , while on the left N ■ 
(1 — Q) nodes the occurring probabilities are 0. 

The secondary-level complexity in the localization 
properties is the deviation of the PDF from the step- 
like function. This deviation corresponds to the shape of 
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the PDF, which can be extracted by using the structural 
entropy [31| . 



N 



Pmlnpm - ln(Q • N). 



(4) 



For the simple step-like condition, we have S s tr = 0. 
S s tr 7^ tells us the shape deviation of PDF from the 
simple one. 

The pair of localization quantities, (Q,S str ), is widely 
used up to date to describe the localization in disor- 
dered and aperiodic systems, and the statistical analysis 
of spectra in diverse fields such as quantum chemistry, 
condensed matter physics, and quantum chaos [32l |33|. 



B. Statistical Properties of the Spectra 

The localization property can also be described with 
the random matrix theory (RMT) [H H [35|. RMT 
is initially developed to understand the energy levels of 
complex nuclei and other kinds of complex quantum sys- 
tems. Recently, the RMT theory has been proposed to 
capture the structure and dynamical properties of com- 
plex networks [ill ]. 

One of the most important quantities in the theory is 
the PDF for the nearest neighbor level spacings (NNLS) 
of the spectrum. It is theoretically and numerically con- 
firmed that at the localization and the extended states 
the PDFs of the NNLS should be Poisson and Wigner- 
Dyson distribution, respectively [37i - [39| . Generally, for 
an intermediate state, the PDF obeys the Brody distri- 
bution, 



U(s) = ^ 



exp 



(5) 



where s is the NNLS and £ the characteristic distribution 
width. The Poisson and the Wigner-Dyson distributions 
are the two extremes with P = 1 and j3 = 2, respectively. 

Introducing the accumulated function, C(s) = 
J U (x)dx, some trivial calculations lead to, 



ln.R(s) 



In 



In 



1 



1-C(s) 



= pins-p\n£. (6) 



From this formula we can testfy the Brody distribution 
and determine reliably the values of the parameters P 
and £. 

To make the spacings s in units of local mean level 
spacing, we should conduct a standard procedure, called 
unfolding. Denoting the spectrum of a network with 
Ai, A2, • • • , Ajv, the accumulation density function for the 
spectrum is G(X m ) = m,m = 1,2, ••• ,N. Fitting this 
relation with a polynomial function, we can separate it 
into a smooth part and a fluctuation part as, 



G(X m ) — G av (X m ) + Gf(X m ). 



(7) 



The NNLS can be obtained as, Sj = G av (Xi + i) — 
Gav(Xi), i = 1, 2, • • • , N — 1. For a complex network, we 
generally have not enough knowledge on its spectrum, 
when the polynomial function fitting method can lead to 
a reliable result. In this paper, the order of the polyno- 
mial function is 17. 



C. Wavelet Transform 

The detailed properties for the PDF of the occur- 
ring probabilities can be used as the measure of the 
global structure symmetries. However, determining this 
PDF is a nontrivial task [4(|. Assume the probabil- 
ity values have been sorted in ascending order, namely, 
p = {p\ < p-2 ■ ■ ■ < Pn}, which can be regarded as the 
profile of the nearest spacing series, Ap = {p2 — pi , pz — 
P2, ■ ■ ■ 1 Pn — Pn-i}- The local structures of Ap can tell 
us the probability distribution function of p. It is found 
that the series p generally behaves multi-fractal. 

The wavelet transform (WT) [4l| is used to detect 
the fractal properties embedded in the ascend-ranked se- 
ries p. The increasing trend in the series p makes the 
box-counting-based techniques invalid to quantify the lo- 
cal scalings. In the wavelet transform, the contributions 
of the polynomial trends can be removed effectively. A 
multi-fractal series can be decomposed into many subsets 
characterized by different local Hurst exponent h, which 
quantifies the local singular behavior and thus relate to 
the local scaling of the series. The statistical properties 
of these subsets can be quantified by the fractal dimen- 
sion D(h) of the subset whose local Hurst exponent is 
h. 

As a standard procedure, we first find the WT maximal 
values, {T g (a, pk(a)) \k = ki,k2, ■ ■ ■ kj}, where a is the 
given scale. The partition function should scales in the 
limit of small scales as, 



Z(a,q) = \ T g( a ' Pk(a))\ 9 ~ a, T 



(9) 



(8) 



k—k± 



The fractal dimension D (h) can be obtained through the 
Legendre transform, 



D(h) = qh - r(q), 



dq 



(9) 



For a mono-fractal structure we have a linear relation, 
T (q) = qH — 1. H is the global Hurst exponent. For 
positive and negative q, r(q) reflects the scaling of the 
large fluctuations and small fluctuations, respectively. 

We use the real analytic wavelet g( n > among the class of 
derivatives of the Gaussian function, by which the poly- 
nomial trends up to n order can be removed. The results 
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FIG. 1: (Color online) The localization quantities (Q, S 8 t r ) 
for the WSSW and BASF networks. There exist complex 
relations between p r or w and (Q, S str ) for the BASF and 
WSSW networks. For the WSSW networks, from p r = to 
p r = 0.02 there exists an abrupt decrease/increase in value 
of Ql S str , as shown in (a)-(b)) respectively. Then with 
the increase of the rewiring probability p r the participation 
ratio tends to increase while the structural entropy tends 
to decrease; (c)-(d) Results for the BASF networks. 



with n = 7 are presented, n = 5 and n — 6 lead almost 
the same results. As comparison, we detect also the scal- 
ing behaviors in the randomized series pr, called shuffled 
series. 

In this paper, we are interested in the character- 
istic point at which the fractal dimension reaches its 
maximum value, (h c , D(h c )). It can tell us the non- 
homogeneous distribution of the series p and the fractal 
characteristics of the principal subset. 



IV. NUMERICAL RESULTS 

We examine the localization behaviors for the cellu- 
lar networks [42j, which are compiled by using a graph- 
theoretical representation of all the biochemical path- 
ways based upon the WIT integrated-pathway genome 
database of 43 species from Archaea, Bacteria and Eu- 
karya (43|. The whole cellular networks consider the cel- 
lular functions as intermediate metabolism and bioen- 
ergetics, information pathways, electron transport, and 
transmembrane transport. The direct edges are replaced 
simply with non-directed edges. We consider only the 
cellular networks with the sizes larger than 500. 

We study also the localization behaviors for the 
the Watts- Strogatz small-world (WSSW) Q and the 
Barabasi- Albert scale-free (BASF) Q networks. For the 
WSSW model, we construct firstly a regular circle lattice 
with each node connecting with its d right-handed near- 
est neighbors. For each edge we rewire it with probability 
p r to another randomly selected node. Self- and double- 
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FIG. 2: (Color online) The relations of S s t r versus Q for 
the WSSW, BASF and whole cellular networks. The lo- 
calization quantities for the distributions, P(p) ~ p-^+v) 
and P{p) ~ -, namely, the critical and localized sates, are 
shown as references. Starting from p(r) ~ r~ a , the set of 
normalized values ofp(-S°),n=l,2,---,A r can be regarded 
as the critical state. Assigning a = 1 ~ 10, the correspond- 
ing values of r] are 0.5 ~ 0.05. N is the size of the considered 
networks. The same procedure is also used to generate the 
localized states by starting from p(r) ~ exp(-pr). The 
localized states with p = 0.01 ~ 100 are generated. The 
localization properties of the BASF networks can be cap- 
tured by the critical state with extremely significantly small 
values of r\. The WSSW and whole cellular networks are 
generally in between the two typical states. 



edges are forbidden. By this way, we can introduce ran- 
domness into the resulting networks. Moreover, com- 
pared with that for the initial regular lattice, the rewiring 
procedure may introduce also "long-range" edges to the 
resulting networks, which can reduce significantly the av- 
erage number of hops required for the nodes to reach each 
other. This is the so-called small- world effect. 

The BASF networks are the results of a preferential 
growth mechanism, which exists widely in diverse fields. 
Starting from several connected nodes as a seed, at each 
growth step a new node is added and w edges are estab- 
lished between this node and the existing network. The 
probability for an existing node being connected with 
the new node is proportional to its degree. Self- and 
double-edges are forbidden. For the resulting networks, 
the number of edges per node obeys a power-law, namely, 
no characterized scale exists in this distribution. 

Fi gure 1 presents the localization quantities, (Q, S s t r ), 
for the networks. For the WSSW networks, the random- 
ness introduced by the rewiring procedure has two com- 
petitive effects, the long-range edges which favors the 
extension, and the broken of symmetry which induces 
the localization. For the BASF networks, the increase 
of w increases the heterogeneity and the connections be- 
tween the nodes, which induce the localization and the 
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FIG. 3: (Color online) The value of Brody parameter j3 
versus network parameters p r and w. (a) WSSW networks; 
(b) BASF networks. 
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FIG. 4: (Color online) The multi-fractal scaling charac- 
teristics of the ascend-ranked series p for the real world, 
the WSSW and the BASF networks. The multi-fractal be- 
havior for the WSSW network with d — 2,N — 2000 and 
p r = 0.05 is presented as a typical example. In the whole 
range of q = — 5 ~ 5, there is only one characteristic point, 
(h c ,D(h c )) = (0.63,0.958). 



extension, respectively. Hence, there exist complex rela- 
tions between p r or w and (Q, S s t r ) for the two kinds of 
networks, as shown in Fig.l(a)-(b) and (c)-(d), respec- 
tively. For the WSSW networks, the participation ratio 
decreases rapidly from 1 to 0.22 whenp r changes slightly 
from to 0.02, and then goes up gradually with the in- 
crease of p r . As for the structural entroy, it increases 
abruptly when p r changes from to 0.02, after that it 
decreases gradually with the increase of p r . 

Figure 2 shows S s tr versus Q. As references, we calcu- 
late also the localization quantities for the critical and the 
localized states, namely, P(p) ~ p^'- 1 ^ and P(p) ~ i 
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FIG. 5: (Color online) The multi-fractal scaling charac- 
teristics of the ascend-ranked series p for the real world, 
the WSSW and the BASF networks. The multi-fractal be- 
havior for the BASF network with w = 8 and N = 2000 
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of q — —5 ~ 5, there is only one characteristic point, 
(h c ,D(h c )) = (0.61,0.942). 
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FIG. 6: (Color online) The branched multi-fractal scal- 
ing characteristics of the ascend-ranked series p for the real 
world and modeling networks. The branched multi-fractal 
behavior for the whole cellular network of M.jannaschii 
is presented as a typical example. The two branches 
q < and q > lead to different characteristic 
points, {h cl ,D{h cl )) = (0.63,0.96) and (h c2 ,D{h c2 )) = 
(0.83,1.03). 



respectively. Starting from p(r) ~ r _cr , we calculate the 
values of p(^),n = 1,2, ••• , N. The resulting normal- 
ized values can be regarded as the localized state. The 
critical sates with a = 1 ~ 10 are calculated, and the 
corresponding values of rj are 0.5 ~ 0.05, respectively. N 
is the size of the considered networks. The same pro- 
cedure can be used to generate the localized states by 
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starting from p(r) ~ exp(— fir) . The localized states with 
/i = 0.01 ~ 100 are generated. 

The localization properties of the BASF networks can 
be captured by the critical states with extremely small 
values of -q. The WSSW and whole cellular networks 
are generally inbetween the two typical (localized and 
extended) states . 

We can find that the PDFs of the NNLS for all the 
networks can be described very well by using the Brody 
distribution in a unified way. The results for the parame- 
ter (3 are shown in Fig. 3. For the WSSW networks, with 
the increase of the rewiring probability p r , the parame- 
ter (3 increases rapidly from 1.02 ± 0.053 at p r = to 
1.95 ± 0.065 at p r = 0.14. For p r > 0.14, /3 are almost 
same, namely ~ 2.0. That is, the representative eigen- 
vector changes from a nearly localized state (p r = 0) to 
an extended state in this interval of p r . For the networks 
with p r > 0.14, the representative eigenvectors are al- 
most perfectly extended. While for the BASF networks, 
with the increase of w, the more edges can induce the sig- 
nificant extensions of the representative states. (3 reaches 
its asymptotic value ~ 1.90. Due to the heterogeneity, 
the BASF networks can not reach a perfectly extended 
state. 

In a considerable wide range of q, the partition func- 
tions behave scale- invariant as in Eq.(8). There are 
three kinds of typical WT transform results. Here we 
present several typical examples. In the whole range of 
q = — 5 ~ 5, the WSSW network with p r — 0.05 and the 
BASF network with w — 8 are multi-fractal with only 
one characteristic point (h c , D(h c ), as shown in Fig. 4 and 
Fig. 5, respectively. Sometimes, the multi-fractal degener- 
ates to mono-fractal. Fig. 6 gives another condition where 
the fractal behaviors can be separated into two branches, 
namely, q < and q > 0. The characteristic points for 
these two branches are not same. That is, the principal 
subsets for the large fluctuations and the small fluctu- 
ations are different. These three conditions are called 
mono-fractal, multi-fractal and branched multi-fractal, 
respectively. 

The scaling properties for the real world networks and 
the modeling networks are listed in Table I. For the mono- 
and multi-fractals, we present simply the global Hurst 
exponent H and the characteristic point (h c , D(h c )), re- 
spectively. For the branched multi-fractal we give the 
scaling characteristics for the two branches q < and q > 
0, which are separated by the division symbol "/". To 
cite an example, for the cellular network M.jannaschii, 
the characteristic point for the branch q < is (0.63, 0.96) 
and that for the branch q > is (0.83, 1.03). It is denoted 
with (0.63, 0.96)/(0.83, 1.03). The results for the corre- 
sponding shuffled series are presented also. We discard 
the networks that the scaling behaviors of the original p 
and the randomized series pu are undistinguishable. The 
sizes of the WSSW and BASF networks are N = 2000. 
N = 1000,3000 and 4000 lead almost same results (not 
shown in Table. 1). 

The WSSW and BASF networks are almost all mono- 



or multi-fractals with the values of h c mainly in the range 
of 0.66±0.05. However, most of the considered real world 
networks behave branched multi-fractal. The Hurst ex- 
ponents lager than 1 and near correspond to non- 
singularity and white noises, respectively. Discarding 
these trivial conditions, we find that the multi-fractal be- 
haviors are embedded in the branches of q > 0. And the 
values of h c are basically in the range of 0.8 ± 0.05. The 
larger values of h c for the large fluctuations in the real 
world networks show us the much more non-homogeneous 
structures of the PDF of p. That is, the PDF of p for the 
real world networks tend to form much sharper peaks at 
different scales. 

It should be noted that in the present paper the 
structure- induced localization is used as a probe of struc- 
ture properties of complex networks. We detect the local- 
ization properties for the WSSW, BASF model networks 
and the cellular networks, but it dose not imply and re- 
quire any localization-related dynamical process (such as 
waves) occurring on the real-world systems. 



V. CONCLUSIONS 

In summary, the probability distribution function of 
the representative eigenvector is proposed to describe the 
localization properties of complex networks. The local- 
ization quantities (Q,S s t r ), the PDFs of the NNLS and 
the wavelet transform are used to capture the characteris- 
tics of the representative state. The nontrivial structures 
of the networks can induce the localizations of the repre- 
sentative states. At the same time, because of the global 
symmetries of the networks, the representative state have 
nontrivial structures rather than the step-like distribu- 
tion. 

The localization quantities (Q, S s t r ) and Brody distri- 
bution parameter /3 can describe the nontrivial localiza- 
tion properties in a global way. The (Q, S s t r ) values tell 
us that the BASF networks with w — 2 are significantly 
localized compared with the WSSW networks with d = 2. 
It is consistent with the conclusions drawn from the re- 
sults of (3. The whole cellular networks have localization 
properties much closer with the WSSW networks. 

The wavelet transform can tell us the details on 
the nontrivial structures of the representative eigenvec- 
tors. The ascend-ranked series p for the WSSW and 
BASF modeling networks and the whole cellular net- 
works behave mono-fractal, multi-fractal or branched 
multi-fractal. The PDF of p tends to form sharp peaks 
at different scales in a self-similar way. 

This kind of property can shed light on the global sym- 
metries due to the general rules in the construction of the 
networks. Hence, it can be employed as a global mea- 
sure of the network structures. Moreover, the structure- 
induced localizations on networks maybe helpful to un- 
derstand the electronic conduction and heat transport 
properties 44] of nanonet materials. 

A closely relevant topic is the diffusion on complex net- 
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works. Kim et al. 45| report for the first time their works 
on quantum and classical diffusion on WSSW networks. 
The Hamiltonian is same as that in the present paper, 
namely e„ = and t mn = 1 in Eq.(f). An electron is 
localized at a randomly selected node at beginning, then 
the diffusion process is obtained by solving the time- 
dependent Schrodinger equation. It is found that the 
" long-range" edges can fasten the diffusion speeds signif- 
icantly, especially at the transition point from p r — to 
p r 0. This is qualitatively in consistent with our find- 
ings of the significant changes of the participation ratio 
and the structural entropy, (Q,S a tr)i when p r increase 
from p r — to p r = 0.02. 

As for the classical diffusion on networks, a very recent 
work reports the first-passage times (FPT) of random 
walkers in complex scale- invariant media fisl |46| . Many 
real-world networks have self-similar structures, and dif- 
fusion on networks can be regarded to a certain degree 
as the diffusion on fractal media, which has attracted 
intensive attentions for its importance in theories and 
potential use in diverse research fields [47j . 



However, we shoud point out that, it is not trivial to 
compare our results quantitatively with these evolution 
processes. Actually our results are obtained from the 
eigenstates in energy representation, while for the quan- 
tum diffusion the initial state of localizing at a randomly 
selected node is a wave packet and the final state should 
be a superposition of the eigenstates in energy represen- 
tation. How to relate the localization with the classical 
diffusion is definitely interesting but not a trivial task. 
Obviously, detailed works on diffusion on complex net- 
works are required to understand the relation between 
localization and diffusion on networks. 
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